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ABSTRACT 

A simple optical to mid-IR color selection, R — [24] > 14, i.e. f v (2A fun)/ f v {R) > 1000, identifies 
highly dust obscured galaxies (DOGs) with typical redshifts of z ~ 2 ±0.5. Extreme mid-IR luminosi- 
ties (L/_r > 10 12-14 ) suggest that DOGs are powered by a combination of AGN and star formation, 
possibly driven by mergers. In an effort to compare their photometric properties with their rest frame 
optical morphologies, we obtained high spatial resolution (0.05 — 0.1") Keck Adaptive Optics (AO) 
if '-band images of 15 DOGs. The images reveal a wide range of morphologies, including: small expo- 
nential disks (8 of 15), small ellipticals (4 of 15), and unresolved sources (2 of 15). One particularly 
diffuse source could not be classified because of low signal to noise ratio. We find a statistically signif- 
icant correlation between galaxy concentration and mid-IR luminosity, with the most luminous DOGs 
exhibiting higher concentration and smaller physical size. DOGs with high concentration also tend to 
have spectral energy distributions (SEDs) suggestive of AGN activity. Thus central AGN light may 
be biasing the morphologies of the more luminous DOGs to higher concentration. Conversely, more 
diffuse DOGs tend to show an SED shape suggestive of star formation. Two of fifteen in the sample 
show multiple resolved components with separations of ~ 1 kpc, circumstantial evidence for ongoing 
mergers. 
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1. INTRODUCTION 

Recent Spitzer Space Telescope 24 fim images of ex- 
tragalactic survey fields hav e revealed extremely dust ob- 
scure d galaxies (DOGs; e.g. iHouck et al.l l2005: Dc v et all 
I2008L iFiore et al.l 120081 ). Defined by very red optical 
to infrared (IR) colors, R - [24] > 14, i.e. /„(24 
Hm)/f v {R) > 1000, DOGs are redder than the typi- 
cal low redshift ultra-luminous in frared galaxy (ULIRG; 
Lir = 10 12 " 13 ; iDev et al.ll2008h . The ~ 9 square de- 
gree Boot es field of the NOAO Deep Wide-Field Survey 
(NDWFS; iJannuzi fc D^[l999h contains ~ 2600 DOGs 
(out of ~ 20000 24 /im sources with a limiting flux den- 
sity of 0.3 mJy). Spectroscopic redshifts of 8 6 DOGs in 
Bootes reveal a < z > = 1.99 with a z — 0.5 ([Dev et al.l 
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2008). At these distances, the implied total IR luminosi- 
ties of the DOGs are Lib > 10 12 ~ 14 , typically in excess 
of low redshift ULIRGS. The extreme luminosities and 
colors suggest rapid AGN accretion and/or intense star- 
formation activity, heavily obscured by dust at rest-frame 
optical and UV wavelengths. 

Optical- IR spectral-energy distributions (SEDs) of the 
DOGs suggest two classes. So called "power-law" sources 
show a continuous rise in flux density to longer wave- 
lengths. Spitzer IRS mid-IR spectra of power-law DOGs 
show s trong silicate absorption, an indicator of AGN ac- 
tivity (jHouck et al.ll2005t IWeedman et al.ll2006[) . In con- 
trast, so called "bump" sources exhibit a rest-frame 1.6 
fim flux excess thought to be produced by the stellar pho- 
tospheres of cooler stars. Spitzer IRS spectra of bump 
sources show st rong PAH features suggestive of rapid star 
formation (e.g. iFarrah et al.ll2008l ; Desai et al. 2009, in 
preparation). At these redshifts, the "bump" appears in 
the Spitzer IRAC bands (3.5 — 8.0 /Ltm). Based on: (1) 
space densities (Dey et al. 2008); (2) clustering strength 
(Brodwin et al. 2008); and (3) mid/far IR SEDs (Pope 
et al. 2008), DOGS, may be linked to the z = 2 sub-mm 
galaxies (SMGs), which are thought to be experiencing 
merger driven star formation rates of up to several thou- 
sand sol ar masses per year . As with ULIRGs at low 
redshift (| Armus et alJl2007f ). most DOGs are likely to 
have some combination of AGN and star formation activ- 
ity. For example, a stu dy of the ave r age (s tacked) x-ray 
properties of DOGs by IFiore et ail (|2008T l showed that 
even lower luminosity DOG s exhibit hard x-r ay sources, 
suggestive of AGN activity. iPope et all (|2008f) . however, 
suggests that fainter DOGs are primarily powered by star 
formation even at x-ray wavelengths. 

While significant progress has been made in deter- 
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mining the bulk properties of the DOGs, less is known 
about the specific triggers of the AGN and star forma- 
tion. High spatial resolution imaging may help iden- 
tify potential triggers. For instan ce, by analogy with 
local ULIRGs ([Sanders et al. 1988) and as might be ex- 
pected from hierarchic al structure formation models (e.g. 
iSomerville et~aT1 l200lh . DOGs may be the product of 
gas rich mergers. High resolution imaging may reveal 
multi-component systems and tidal features typical of 
mergers. Alternatively DOGs may be the product of the 
monolithic collapse of primordial clouds s een at the time 
of first assembly (e.g. lEggen et al.l [T962L If the DOGs 
represent the formation of today's most massive systems 
(typically elliptical galaxies), they may alrea dy preferen- 
tially exhibit elliptical morphologies at z = 2 fZirm et al.l 
[2007T). 

iMelbourne et alj (|2008bl hereafter Mel08) and Buss- 
man et al. (2009, in preparation) describe our early 
efforts to obtain high spatial resolution Keck Adaptive 
Optics, (if' -band) and Hubble Space Telescope (HST; 
NICMOS iJ-band; ACS i-band; and WFPC2 U-band) 
images of the DOGs. These studies focussed on power- 
law sources and, especially at longer wavelengths (H and 
if '-bands), typically found smooth systems with small 
sizes (i.e. typically R\/2 < 3 kpc). The HST optical 
data (rest UV) tended to contain more substructures as 
might be expected from multiple star forming regions. 

Unfortunately, there were significant limitations to 
these first efforts. The Keck Adaptive Optics (AO) imag- 
ing in particular was a very small sample size, contain- 
ing only three DOGs. The AO sample were among the 
brightest in the entire NDWFS sample, none with a 
known redshift. While the sample size for the HST study 
was significantly larger, 30 systems, the HST data also 
had limitations. The very faint optical (rest UV) fluxes 
for the DOGs resulted in low signal-to-noise ratio (S/N) 
in the optical HST imaging, especially for the WFPC2 
data. The HST NICMOS (TJ-band, rest blue) images of 
the DOGs were higher S/N than the optical data, but the 
NICMOS ii-band imaging had significantly lower spatial 
resolution (~ 0.15") than the HST optical or Keck AO 
if-band (~ 0.05") data. 

This follow-up paper presents a significantly larger 
sample of DOGs with Keck AO imaging in the if -band, 
now 15, including 2 with overlap with the HST sample. 
The sample now spans the range of 24 /im flux densi- 
ties for NDWFS DOGs ( 0.3 mJy < f„(24 fxm) < 5 
mJy), and includes 6 potential bump sources. Section 
3 gives the Keck AO derived morphologies of these sys- 
tems and correlates them with other observed properties 
of the sample. Results are discussed in Section 4. 

Throughout, we report Vega magnitudes and assume 
a A cold-dark-matter cosmology: a flat universe with 
H = 70 km/s/Mpc, Q m = 0.3, and fl A = 0.7. 

2. THE DATA 

2.1. Sample Selection 

Adaptive optics makes use of a reference source to 
track and correct for atmospheric turbulence. AO cor- 
rected images approach the diffraction limited resolu- 
tion of the telescope which for Keck in the if -band is 
~ 0.05", similar the the resolution of HST in the op- 
tical. Natural Guide Star (NGS) facilities use a star 
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Fig. 1. — Optical/IR color-magnitude diagram for DOGs in 
Bootes (points). The DOGs in the AO sample are shown as dia- 
monds. The data show both the sharp 24 fim flux limit and the 
color definition for DOGs. The AO sample spans the full range of 
DOGs in color-magnitude space. 

near the science target as a reference source. NGS guide 
stars need to be fairly brigh t (typically R < 14.5 for 
Keck, IWizinowich et all 120001 : Ivan Dam et~all I2004D for 
the AO system to work. In addition, the stars need to 
be near (typically within 30") the science target because 
the AO correction falls off with distance from the guide 
star (anisoplanetism). A Laser Guide Star (LGS) facil- 
ity uses a laser to produce a reference source high in 
the Earth's atmosphere for the AO correction. The laser 
spot is used to track the high order wavefront errors pro- 
duced by the turbulent atmosphere. LGS systems still 
require a "tip-tilt" reference star to correct the low order 
(tip/tilt) terms and eliminate image motion. However 
these tip-tilt stars can be fainter and significantly further 
from the science targe t than NGS guide stars (R < 17.5 
and d < 55" for Keck. IWizinowich et al.l l2006). Another 
benefit of the LGS facility is that the laser spot is at 
the location of the science target, reducing the affects of 
anisoplanetism. Of the ~ 2600 DOGs in the Bootes field 
of the NDWFS ([Jannuzi fc Devi[l99l . roughly a third 
(918) are observable with with Keck LGS AO (using tip- 
tilt reference stars with R < 17.5 [mag], at an angular 
separation d < 55"). Only 15 are observable with Keck 
NGS AO (using NGS guide stars with R < 14.5 [mag], 
at an angular separation d < 30"). 

A wide variety of selection criteria were used to select 
the DOGs in this paper. Fo urteen were chosen fro m the 
Bootes field of the NDWFS ([Jannuzi fc Devi[l999h . One 
was chosen from the Spitzer First Look Survey (FLS; 
lYan et al.l [2007h . In 2007, the three LGS targets were 
selected to provide the best possible S/N and AO perfor- 
mance, i.e. among the brightest DOGs in the if -band, 
and close to a bright AO ti p-tilt star. Initia l result s 
on these sources were given in IMelbourne et all |2008b). 
Three NGS targets were selected to be near to bright 
NGS guide stars, with no selection for if -band bright- 
ness. None of the DOGs observed in 2007 had spectro- 
scopic redshifts. 

The 2008 AO selection criteria were relaxed to primar- 
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Fig. 2. — The 2007 Keck AO observations of 6 DOGs (left) and their associated PSFs (middle). Contours are over-plotted to demonstrate 
the differences between the galaxy and PSF. ID radial intensity profiles of the DOGs (solid line) and PSFs (dotted line) are shown (right). 
The DOGs tend to have more extended profiles than their associated PSFs. The box size is ~ 1.5" on a side. DOGs 1-3 were observed 
with the narrow-field camera (pixscale = 0.01" pix" 1 ) in LGS mode, while DOGs 4 - 6 were observed with the wide-field camera (pixscale 
= 0.04" pix" 1 ) in NGS mode. 



4 




Fig. 3. — Same as Figure [2] only now for 2008 observations. DOGs 7-9 were observed with the wide field camera (pixscale = 0.04" pix ) 
in NGS mode, while DOGs 10 and 11 were observed in the narrow-field camera (pixscale = 0.01" pix -1 ) in LGS mode. 
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Fig. 4. — Additional 2008 observations made with the wide field camera (pixscalc 
mode, while DOG 15 was made in the NGS mode. 

ily select DOGs with existing spectroscopic re dshifts, in- 
cluding two NGS targ ets and five LGS targets (|Yan et al.l 
l2007tlDev etalll2008h . Two additional NGS targets ob- 
served in 2008 did not have spectroscopic redshifts. Com- 
plicating the 2008 selection criteria, two DOGs were se- 
lected because of interesting morphology in existing op- 
tical and near-IR HST imaging. 

While the selection criteria are non-uniform, the final 
sample of 15 DOGs spans the full range of 24 /xm fluxes 
and R — [24] colors of the larger Bootes DOG sample 
(Figure Q}. 

2.2. AO Observations 

Keck AO observations were obtained for 15 DOGs over 
the course of 5 nights in the Spring of 2007 and 2008. 
Table Q] details the observing conditions on these dates. 
Because the first two nights in the 2007 campaign had 
heavy extinction by clouds (over 1 magnitude of extinc- 
tion in R) 1 we could not propagate the laser. Instead we 
observed 3 NGS targets. Conditions were better during 
the final 2007 observing night, which produced images 
of 3 DOGs from the LGS target list. In 2008, we ob- 



: 0.04" pix" 1 ). DOGs 12 - 14 were made in the LGS 



TABLE 1 

Keck Observation Run Summary 



Obsrun 


UT Date 


Run Type 


Conditions 


1 


May 12, 2007 


NGS 


thick cirrus 


2 


May 21, 2007 


LGS 


scattered clouds 


3 


May 19, 2008 


NGS 


clear 


1 


May 20, 2008 


LGS 


clear 


5 


June 2, 2008 


LGS 


clear 



served 4 DOGS with the NGS facility and 5 with the 
LGS facility, all under good observing conditions. For a 
more detailed description of our observing strategies w ith 
the Keck LGS AO facility see lMelbourne et all (|2008aU rJ). 

Observations were made with the NIRC2 infrared cam- 
era in the K' filter. Details of the observations including 
exposure times, pixel scale, and estimated AO perfor- 
mance (full- width-half- max of the point-spread-function) 
are given in Table [2] The brightest sources were observed 
with the narrow (10" x 10") camera, which has the finest 
pixel sampling (0.01" pix -1 ), taking full advantage of the 
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TABLE 2 

Observing Summary 



















per pix 


DOT? 

For 


TipTilt 


or INCjo 


object 


/~~\ rr; : 1 ~YT a 

Omcial JNarae 


z 


f24 


R - [24] 


Obsrun 


exptimc 


pixscale 


bis limit 


1 WHM 


R 


scp 








[mj] 


[mag] 




[min] 


["/pix] 


[mag/(") 2 ] 


["] 


[mag] 


["] 


DOG 1 


SST24 T143234 94-333637 




2.924 


14.47 


2 (LGS1 


30 


0.01 


17.4 


0.057 


15.0 


26 


DOG 2 


SST24 1142801 (1+341525 




2.487 


15.44 


2 (LOSI 


30 


0.01 


17.4 


0.053 


14.3 


20 


DOG 3 


SST24 11 42944 94-324332 




1.148 


14.09 


2 (LGS^ 

1 UVJU 1 


9 


0.01 


16.9 


0.053 


14.5 


27 


DOG 4 


SST24 J143117.1+332024 


_ 


0.344 


14.64 


1 (NGS) 


60 


0.04 


18.5 


0.088 


11.6 


28 


DOG 5 


SST24 J142825.5+343830 




0.310 


15.52 


1 (NGS) 


68 


0.04 


19.3 


0.088 


13.7 


22 


DOG 6 


SST24 J142616.2+352116 




0.337 


14.30 


1 (NGS) 


60 


0.04 


19.0 


0.088 


13.0 


29 


DOG 7 


SST24 J142925.9+345151 




0.359 


14.45 


3 (NGS) 


96 


0.04 


19.4 


0.104 


14.2 


23 


DOG 8 


SST24 J143032.8+340046 




1.854 


15.56 


3 (NGS) 


81 


0.04 


19.1 


0.072 


14.0 


23 


DOG 9 


SST24 J143641. 1+350207 


1.948 


0.332 


14.64 


3 (NGS) 


63 


0.04 


19.2 


0.072 


14.2 


24 


DOG 10 


SST24 J143335.6+354243 


1.297 


5.577 


14.06 


4 (LGS) 


24 


0.01 


17.2 


0.052 


15.5 


47 


DOG 11 


SST24 J143027.2+344008 


1.370 


1.169 


15.28 


4 (LGS) 


39 


0.01 


17.4 


0.051 


16.8 


41 


DOG 12 


SST24 J143025.7+342957 


2.545 


2.471 


15.23 


4 (LGS) 


69 


0.04 


19.3 


0.104 


14.9 


30 


DOG 13 


SST24 J142538.1+351856 


2.260 


0.846 


17.60 


5 (LGS) 


116 


0.04 


19.5 


0.070 


14.3 


19 


DOG 14 


SST24 J143424.5+334542 


2.263 


0.861 


15.56 


5 (LGS) 


36 


0.04 


19.0 


0.078 


13.8 


34 


DOG 15 


MIPS 16113 


1.930 


1.042 


14.20 


3 (NGS) 


66 


0.04 


19.1 


0.075 


12.9 


31 



a iHouck eFUl |2005|) 

AO correction. Fainter DOGs, not detected in sky sub- 
tracted 5 minute exposures with the narrow camera, were 
observed with the wide field camera (0.04" pix -1 ) which 
has a field of view large enough (40" x 40") to contain 
one or more bright sources to align successive images for 
stacking. 

Individual images were typically between 2 and 5 min- 
utes. A telescope dither was applied between successive 
images. Total exposure times varied from ~ 30 minutes 
for the brightest sources to over an hour for the faintest 
sources. One source, DOG 3, was only observed for 9 
minutes because of an instrument fault. Final reduced 
images were created from a clipped mean of the individ- 
ual exposures after subtracting sky background, dividing 
Hatfield variations, and correcting for camera distortions. 
Sky and flatfield frames were created from the actual sci- 
ence images with sources masked out. Table[2]records the 
per pixel surface brightness limit of each image. Note, 
images taken with the narrow camera have brighter per 
pixel surface brightness limits, but have 16 pixels for each 
wide camera pixel. 

AO images of the DOGs are shown in FiguresOGS and 
[4j Contours are over-laid to show the distribution of flux 
more clearly. To determine if the DOGs are resolved, 
the flux profiles of the DOGs are compared to AO point- 
spread-functions (PSFs) observed nearby in time to the 
target observations. The PSFs are discussed below. 

2.3. AO Point- Spread- Function 

The AO PSF is a product of the atmosphere and the 
corrective optics, both of which are changing on short 
time scales. Not only does the AO PSF change tempo- 
rally it also changes spatially across the field. This is 
because the AO correction drops with separation from 
the both the tip-tilt star (anisokinetisism) and the laser 
spot (anisoplanatism). Accurate understanding of the 
AO PSF is necessary for the morphological analysis that 
follows. Therefore, in addition to observing the science 
targets, we also attempted to image stars, representa- 
tions of the real-time AO PSF. 

For DOGs 1, 4, and 12 these PSF stars were actually 
in the science fields, at a similar separation from the 
AO guide star as the science target. For several of the 



NGS targets (DOGS 7, 8, 15) we observed PSF star pairs 
where the observed star had a similar separation from its 
guide star as the actual science target. For most of the 
LGS targets (DOGs 2, 10, 13, 14) we observed the tip- 
tilt guide star as a measure of the AO PSF. Because of 
anisoplanatism, observing the AO guide star is generally 
not a good measure of the PSF for NGS observations. It 
is a more reasonable estimate for LGS observations where 
anisoplanatism is minimal and the isokinetic angle (the 
angular scale over which tip-tilt varies) is large. For those 
DOGs without direct measure of a PSF (DOGs 3, 5, 6, 9, 
11) we use the PSF of an alternative DOG taken nearest 
in time to the science observation. Because the PSF is 
generally elongated in the direction of the laser spot or 
NGS guide star, we rotate the PSFs to best match the 
expected PSF at the position of the DOG. 

The best estimate of the effective resolution of the 
AO images at the positions of the DOGs are given by 
the the full-width-half-max (FWHM) of the PSFs, pre- 
sented in Table 2. For the wide camera data, the typical 
FWHMs are ~ 0.08". The typical narrow camera data 
have FWHM measures of 0.05 - 0.06". If the DOGs con- 
tain multiple point-like structures with separations larger 
than these limits they will be resolved in our images. 

3. MORPHOLOGY MEASUREMENTS 

FiguresdlGIl andQ]show images of the DOGs and their 
associated PSFs. A comparison of the DOG profiles with 
the PSF profiles suggest that the bulk of the DOGs are 
resolved (13 of 15). DOGs 3 and 11, however, appear to 
be unresolved, or only marginally resolved. 

From these images, we make three types of morpho- 
logical measurements; (1) circular aperture photome- 
try to estimate half light size, and concentration (de- 
fined below); (2) 2D galaxy profile fitting with GALFIT 
(|Peng et al.ll2002l ) to measure a PSF corrected effective 
radius and Sersic index (useful for differentiating between 
disk-like and elliptical-like systems, see Section 3.2 for a 
definition of the Sersic profile); and (3) multiple compo- 
nent fits with GALFIT to reveal any sources that show 
either a strong central point source indicative of AGN 
activity, or sources that show two or more resolved com- 
ponents suggestive of merging. These measurements are 
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discussed in detail below. 

3.1. Circular Aperture Photometry 

For each galaxy we measure circular aperture photom- 
etry about the center of the source, where the center is 
se lected by a Gaussian fit to the galaxy. As was done 
in iBershadv et all (|2000h we use a curve of growth tech- 
nique to determine the total flux in each system. We use 
a maximum aperture size of 2", 5-10 times larger than 
the typical DOG half-light size. From total flux, we de- 
rive the radius which contains half the light, r^o, and the 
galaxy concentration. Concentration is defined as: 

C = 5log (r so/ r 2 o), (1) 

where rgo is the radius that contains 80% of the light, 
and r2o the radius c ontaining 20% of the light (e.g. 
IBershadv et alj 120001 : IConse licc 2003) . Concentration 
has been shown to correlate with galaxy Hubble type. 
Typical concentrations for galaxies in the local universe 
range from C ~ 3 for late ty pe disk galaxies, to C ~ 4 
for early types and ellipticals (iBershadv et al.|[2000f ). 

For the DOGs, typical concentrations range from C ~ 
4 at the high luminosity end to C < 3 at the low luminos- 
ity end. Half-light sizes vary from roughly 0.2 — 0.5" or 
1 to 4 Kpc (assuming the DOGs are at z = 2). Measure- 
ments of DOG half-light radius and concentration are 
given in Table [31 We do not correct C and r$o for PSF 
effects. However, for comparison, we provide measure- 
ments of C and r$o of the PSF stars (Table [3]). Typical 
PSF concentrations are C ~ 3.9, and half-light sizes are 
r 50 ~ 0.15". 

Cloning local galaxi es into the high redshift universe, 
IBershadv et al.l (|2000D showed that C is robust (mea- 
sured to within 10 — 20%) for galaxies with half-light 
radii larger than two resolution elements (i.e. 4 pixels). 
In most cases, the DOGs have half-light sizes larger than 
4 resolution elements, suggesting that their C measures 
are robust. For the bulk of the DOGs in our sample, the 
estimated uncertainty in C due to small half-light sizes is 
significantly larger than that introduced by photometric 
errors. We combine these two uncertainties as a final es- 
timate of the uncertainty on C. The uncertainties in the 
half-light sizes and concentrations range from 10-20%. 

3.2. Single Sersic Profiles 

Given a galaxy image and its associ ated PSF, the 2D 
galaxy profile fitting routine GALFIT (|Peng et al.ll2002f) 
uses a chi square minimization routine to estimate the 
best-fit Sersic profile of the galaxy. The Sersic profile is 
defined by: 

£(r) = Eeezph^-) 1 /" - 1)], (2) 

where S is the galaxy surface brightness, n is the Ser- 
sic index, and r e is the effective radius. All profiles 
are assumed to be axially symmetric ellipses with r = 
(x 2 + y 2 /q) 1 / 2 , and q equal to the minor to major axis 
ratio of the ellipse. A Sersic values of n = 1 reduces 
to an exponential profile typically associated with disk 
galaxies. A Sersic value of n — 4 is a de Vaucouleurs 
profile, typically associated with elliptical galaxies. For 
each DOG, GALFIT fits for total flux, Sersic index, ef- 
fective radius, semi-major to semi-minor axis ratio, and 



position angle. GALFIT convolves each test model with 
the AO PSF, and minimizes the residual difference of the 
model image with the actual data. Because GALFIT has 
trouble disentangling sky from galaxy light for the lower 
surface brightness edges of the galaxies, we provide GAL- 
FIT with an independent measure of the sky. We use the 
median of the pixels in a 4" box surrounding the galaxy, 
with pixels associated with the galaxy removed. 

The typical DOG single Sersic profile has a Sersic index 
of n < 2 (9 of 14 that were measurable), although one 
of these systems, DOG 11, was consistent with a point 
source (i.e. effective radius r e // < 0.05"). Several had a 
Sersic index n > 3 (5 of 14), although 1 of these systems, 
DOG 3, was also consistent with a point source. DOG 6 
was too low surface brightness to measure with GALFIT. 
The model parameters, Sersic index, and effective radius, 
are recorded in Table [3] Because the half-light sizes of 
the DOGs are small, typically only 50% - 100% larger 
than the PSF half-light sizes, the effective radii measured 
by GALFIT (which corrects for the PSF) are typically 
smaller than the DOG half-light sizes calculated with 
aperture photometry and no correction for the PSF. 

Figure [5] shows the DOG images (left), the best fit sin- 
gle Sersic model for each DOG (middle), and the residual 
difference between the two (right). For the majority of 
DOGs, these residuals show minimal structure. How- 
ever, several show significant structure possibly indicat- 
ing a central point source or multi-component system. 
Most notable of these are DOGs 1, 10, and 13. We will 
examine these structures in more detail in the following 
sections. 

While by-eye examination of the residuals is the pri- 
mary method for determining the adequacy of the fits, 
we also calculate a quality statistic, Q, given by: 

stdev of residuals within 0.4" of the galaxy center 
stdev of pixels in the neighboring sky 

(3) 

For a good fit, Q approaches 1. This statistic was 
chosen because the dominant source of uncertainty was 
the infrared thermal sky flux. The sky background can 
change dramatically on short time scales, and it was 
removed during image processing. By comparing the 
GALFIT residual image to noise in the sky, Q more ad- 
equately describes the quality of the fit, than the GAL- 
FIT chi-square measures which do not account for the 
removed sky background. For each GALFIT result, Ta- 
ble [3] provides a Q measure. For the single Sersic fits, 
typical Q values range from Q= 1 to 2. 

3.3. Sersic + PSF Profiles 

We also use GALFIT to estimate contribution of cen- 
tral point sources to the DOG light profiles. An un- 
resolved point source in the center of a galaxy could 
indicate the presence of an AGN, unresolved bulge, or 
unresolved central star burst with a size less than 1 kpc. 
In this mode, GALFIT simultaneously fits for a Sersic + 
PSF profile. Table [3] gives the total flux ratio of the best 
fit central point source to the rest of the Sersic compo- 
nent. We caution that while GALFIT produces a mea- 
surement, these two component models are not necessar- 
ily better than a one component fit. 

To interpret the results, we compare the two compo- 
nent fits to the single Sersic fits, looking for improve- 




Fig. 5. — GALFIT models for the 15 DOGs . The left column shows the actual science data. The middle column shows the best fit single 
Sersic model from GALFIT (Peng ct al. 2002). The right panel shows the residual difference between the two. Images are ~ 1.5" on a side. 



ments in Q and in the residual images. For instance, 
DOG 10 shows a large improvement in Q when adding 
an additional point source component. As will be dis- 
cussed in the next section DOG 10 is actually best fit by 
a 3 component system Sersic + Sersic + PSF (See Figure 
(HI). The point source component is roughly ~ 18% of the 
luminosity of the Sersic component. 

DOG 12 also appears to host a strong central point 
source. For DOG 12 the point source has a PSF/Sersic 
flux ratio of 80%. The single Sersic fit for DOG 12 had a 
Sersic index of 3.47, consistent with an elliptical galaxy. 
However, after including a PSF in the fit, GALFIT finds 
that the extended component has a Sersic index of 0.55, 
more consistent with a disk. This is the only galaxy in 
the sample for which including a point source in the fit, 
fundamentally changes the Sersic index from elliptical- 



like to disk-like. 

Four other DOGs have PSF/Sersic flux ratios greater 
than 20%, (DOGs 3, 4, 8, and 11). However, in the case 
of DOG 4, which has the largest PSF/Sersic flux ratio 
of 2.6, the Q value and the residual image is actually 
worse for the 2 component model than for the 1 compo- 
nent model. Therefore, we do not believe that the two 
component, PSF dominated, model is an accurate repre- 
sentation of the galaxy. 

For DOGs 3 and 11, the single Sersic fits were already 
consistent with point sources. Recall, the single com- 
ponent Sersic fits to DOGs 3 and 11 had effective radii 
of 0.02" and 0.03" respectively, i.e. the size of PSFs. 
Therefore, any additional components were difficult to 
measure. When GALFIT fit for an additional compo- 
nent in the DOG 3 model, it measured a nonsensical size 
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Fig. 6. — Different GALFIT models of DOG 10 (middle column). The science image is shown in the left column and the residual difference 
between the data and the model is shown in the right hand column. The top row uses a single Sersic model, which is not a good fit to the 
galaxy light. The middle row uses the best fit 3 component model: Sersic 1 + PSF + Sersic 2. This model shows minimal residuals. The 
bottom row uses the same 3 component model, only Sersic 2 is not subtracted off. The two resolved components are separated by only 
0.13" (~ 1 kpc at this redshift), suggesting that DOG 10 is undergoing a merger. The box size is ~ 1.5" on a side. 

of 90 arcsec for it (see Table [3]) ; usually this indicates 
that GALFIT has confused sky with object. Therefore 
the single component, point source dominated fit is a 
better approximation for the morphology of DOG 3 than 
the two component model. In the case of DOG 11, the 
GALFIT two component Sersic + PSF model has a Ser- 
sic component with an effective radius the size of a PSF. 
Thus GALFIT effectively chose to fit the galaxy with a 
PSF + PSF model. This again suggests that the single 
Sersic fit is more appropriate than a two component fit. 
Thus, both DOGs 3 and 11 appear to be point-like and 
any underlying resolved galaxy component is too faint to 
detect. 

In summary, five DOGs show strong evidence for lumi- 
nous, central, point-like structures (DOGs 3, 8, 10, 11, 
and 12). GALFIT finds less convincing evidence for lu- 
minous point-like structures in the remaining 9 systems 
that were measurable. 

3.4. Multiple Resolved Components 

In the third round of GALFIT modeling, we exam- 
ine each DOG for an additional resolved component. In 
Mel08, we showed that DOGs 1-3 did not show evidence 
for a second offset resolved component such as would be 



seen in an ongoing merger event. Mel08 also showed ex- 
amples of what a second resolved component would ap- 
pear like in the GALFIT residual images. In this larger 
data set, DOGs 10 and 13 show residual patterns sug- 
gestive of a second resolved component. 

The best fit model for DOG 10 is actually a 3 com- 
ponent model: the main galaxy (Sersic 1), an unresolved 
point source associated with the main galaxy (PSF), and 
a second resolved component (Sersic 2). The total flux 
ratio of Sersic 2 to Sersic 1 is ^20%. If any of these 
components are not included in the fit, there is signifi- 
cant structure in the residuals, and the quality code, Q, 
degrades. Figure [6] shows the difference between a sin- 
gle Sersic fit, and a Sersic 1 + PSF + Sersic 2 fit. Also 
shown is the residual without subtracting off the second 
resolved component. Clearly a second resolved compo- 
nent is contributing to the light profile of DOG 10. 

DOG 13 is also best fit by a 3 component model, Ser- 
sic 1 + PSF + Sersic 2 (Figure 0. In this system, the 
component that contains the bulk of the light is low sur- 
face brightness, and not obvious in the image until the 
higher surface brightness component is subtracted off (to 
best see the low surface brightness component look at the 
residuals in the bottom row of Figure [7]) . It has a large 
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Fig. 7. — Same as Figure l6l only now for DOG 13. The single Sersic model (top row) is not a good fit. The three component model, 
Sersic 1 + PSF + Sersic 2, is the best fit model (middle row). The bottom row shows the same as the middle row, only the main galaxy 
is not subtracted. This third panel is surprising. The main galaxy is low surface brightness and not obvious until after subtracting off the 
companion, and yet it has 5 times the total flux of the companion. The PSF is associated with the main galaxy. Note, this set of images 
is zoomed out from the previous images to show a larger area (3" on side). 
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0.51 


4.27 c 


0.88 


0.50 


1.19 


0.70 


0.53 


3.98E-02 


1.16 










DOG 8 


3.6 


0.13 


3.5 


0.31 


2.60 c 


4.77 


0.21 


1.44 


2.26 


0.30 


2.27E-01 


1.38 










DOG 9 


3.6 


0.13 


2.9 


0.30 


2.55 


1.61 


0.21 


1.04 


0.98 
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0.41 
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a Concentration Sersic Indcx C Assuming z = 2 (stdcv of residuals in 0.4" radius aperture about the center of the galaxy) / (stdev of pixels in the neighboring sky) e (flux of point 



source) / (flux of Sersic) (flux of fainter Sersic) / (flux of brighter Sersic)^ This is a very large size and probably un-physical. This galaxy appears to be point source dominated 
and after subtracting off this point source GALFIT struggles to fit another component.' 1 For the two unresolved DOGs, Rj/2 is an upper limit on the physical size. 
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Fig. 8. — A wide field view of DOG 9 reveals a second object 
with a projected separation of 1.6" (~ 13 kpc if both at z = 2). 

half-light size, ~ 0.49", and contains about 5 times the 
total flux of the higher surface brightness component. 
The central PSF component contains only 2% of the flux 
of the main galaxy. There also appear to be several ad- 
ditional nearby structures not included in the fit - e.g. 
a faint structure appears about ~ 0.8" to the south-east 
of the main galaxy. It is not clear if these are additional 
components in the same system or just background noise. 
Several other DOGs may have companions. 

• Figure [8] reveals that DOG 9 has a neighbor, with 
similar flux, at a separation of 1.6" (~ 13 kpc if 
both at z = 2). We do not know if this companion 
is at the same redshift as DOG 9. 

• DOG 11, which is point source dominated, shows 
residuals suggestive of a very faint companion 1/10 
of the flux of the main system. (See Figure 5. The 
possible companion is a fuzzy patch to the north 
of the main galaxy in Figure 5.) Because the flux 
ratio is so small it is not clear if this is truly an 
additional system. 

• DOG 6 is very diffuse and GALFIT was not able 
to fit the galaxy, but the residuals (Figure 5) after 
subtracting off the sky suggest that it may be made 
up of faint multiple components. 

• DOG 1 also shows significant structure in the GAL- 
FIT residual image (Figure 5). This is best ac- 
counted for by adding a second component to the 
fit. This second component is consistent with an 
exponential bulge, rather than a merging compan- 
ion. 

For the DOGs that show two resolved components Ta- 
ble [3] gives the total flux ratio of those two components, 
and their separations from the main galaxy. 

4. DISCUSSION 

Mel08 analyzed DOGs 1-3, among the most luminous 
in the sample. Compared with other high redshift galaxy 
samples, the 3 DOGs in Mel08 showed high concentra- 
tions and small half-light radii. For example, radii were 
significantly smaller than the sizes of z = 1 luminous 
infrared galaxies (LIRGs, Lib, ~ 10 11-12 Lq) observed 
with Keck AO by th e Center for Adaptive O ptics Trea- 
sury Survey fCATS. [Melbourne et al.ll200 8a). Concen- 
trations were higher than both the L IRGs and z = 2 
sub-mm sources fe.g. lPope et~aHl200l . 

In this work, the significantly larger sample spans the 
full range of 24 micron flux densities of DOGs in Bootes. 
Although it contains a higher percentage of DOGs at 
the high flux end, the sample includes DOGs that are al- 
most order of magnitude fainter than the DOGs in Mel08. 
While many of the brighter DOGs show high concentra- 
tions (similar to Mel08), as evidenced by radial profiles 



only marginally more extended than PSFs (Figure [2]), 
several of the fainter DOGs show substantially more dif- 
fuse morphology. 

4.1. A Morphology - Luminosity - Color Relation 

Figure [5] plots galaxy concentration (at 2.2 /im) vs. 
24 /xm flux density. While flux is an observed quantity, 
DOGs in Bootes have been shown t o have a tight re dshift 
distribution (< z >~ 2, a z ~ 0.5. lDev et alJl2008f ). and 
therefore flux should be closely tied to physical luminos- 
ity. Galaxy concentration appears to be correlated with 
24 /im flux. Brighter DOGs have higher concentrations 
than the fainter ones. A chi squared test, fitting concen- 
tration vs. 24 /im flux density, rules out a constant fit 
to the DOG data (i.e. no correlation) at the 99% level. 
A linear relation between the two parameters can not be 
ruled out by the chi squared test. Also shown in Figure 
1 are z ~ 1 LIRGs from the CATS sample. The LIRGs 
appear to continue the trends seen in the DOG sample. 

There is some concern that the statistically significant 
correlation between galaxy concentration and 24 (im mi- 
cron flux is not physical but rather the result of the way 
the data were taken. The bulk of the low 24 /im DOGs 
were observed in NGS mode which, at large separations 
from the guide star, can result in a significantly poorer 
AO correction compared with LGS data. However, there 
is not a statistically significant correlation between PSF 
concentration and galaxy concentration. Likewise a chi- 
square test of a fit between PSF concentration and galaxy 
24 /im flux cannot rule out a constant fit (i.e. no correla- 
tion). Therefore the correlation between galaxy concen- 
tration and 24 /im flux appears to be real and unrelated 
to deficiencies in the data. 

Figure[H]also shows the half-light sizes (measured at 2.2 
/zm) of the galaxies plotted vs. 24 /im flux density. The 
brightest systems tend towards smaller sizes compared 
with the fainter systems. As was shown in Mel08, the 
brighter sources are smaller than the z ~ 1 LIRGs. The 
fainter DOGs, however, tend to have sizes comparable to 
the LIRGs. 

The spectral energy distributions of the DOGs may 
help explain these trends. Figure [TU] shows concentra- 
tion plotted as a function of rest-frame near-IR color 
from Spitzer IRAC photometry. In these plots, rest- 
frame near-IR color is used as a proxy for SED shape. 
DOGs with power-law SED shapes (i.e. rising flux den- 
sity to longer longer wavelengths) should have log total 
flux ratios less than in both of the these plots. However, 
SEDs with a strong flux excess from cool stellar atmo- 
spheres (i.e. "bump" sources) should show log total flux 
ratios greater than in one or both of these two plots. 
Figure ITD1 shows that the six most diffuse DOGs have a 
potential stellar bump. Four (DOGs 4, 5, 9, and 13) have 
a flux excess at 4.5 microns (left), two of which (DOGS 
4 and 5) also have a flux excess at 5.8 microns (right). 
The two most diffuse DOGs (DOGs 6 and 7) have a flux 
excess at 5.8 microns. The other nine DOGs do not have 
SED shapes suggestive of bump sources. Their SEDs rise 
into the IR as expected for power-law sources. 

As has been discussed in previous papers, DOGs 
with power-la w SEDs are associated with significant 
AGN activity dHouck et al.ll2005t iWeedman et "all 120061 : 
iBrand et alj I2007T K On the other hand, star formation 
dominates the mid-IR flux output of DOGs with rest- 
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Fig. 9. — Concentration (left) and half-light radius (right) as a func tion of 24 fim fl ux density. While 24 (im flux is not an intrinsic 
quantity, DOGs have a tight redshift distribution (< z >~ 2, <j z ~ 0.5, Dey et al. 2008), and therefore 24 fim flux should be closely tied 
to rest-frame mi d-IR luminosity. DOGs a re shown as filled stars. A comparison sample of z = 1 LIRGs from the the CATS survey is also 
shown (triangles; Melbourne et al. 2008a). For the DOGs, Concentration appears correlated with 24 fim flux. A chi square test rules out, 
at the 99% level, a constant fit (i.e. no correlation) to the C vs. 24 fim data. A similar trend is seen for half-light size, with the more 
luminous systems tending to show smaller sizes. The more luminous DOGs are typically smaller with higher concentrations than the 2 = 1 
LIRGs. 
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Fig. 10. — Concentration vs. IR color (from Spitzer IRAC imaging; Eiscnhardt et al. 2004). /„(4.5 fim) / /„(5.8 fim) is shown on the 
left, while /„(5.8 fim) / /„(8.0 fim) is shown on the right. "Bump" sources, with a rest 1.6 fim SED stellar bump are likely to have a 
log color ratio greater than in at least one of these two plots (the six bump candidates are circled). Power-law sources, more likely to 
be powered by AGN, should show log color values less than in both these plots. Interestingly, the most diffuse systems appear to have 
color terms suggestive of SED bumps. The two most diffuse sources appear to show a bump at 5.8 ^tm(right). The potentially 4 additional 
bump sources at 4.5 fim (two of which appear bump like at 5.8 fim as well), also exhibit diffuse morphologies (left). 

mation dominated, are more diffuse. If AGN exist in the 
bump DOGs, they appear to be sufficiently enshrouded 
so as to not bias the rest-frame optical morphologies to 
higher concentrations. Of the six potential bump sources, 
GALFIT finds none with PSF/Sersic ratios greater than 
12%, except for DOG 4, where we suspect that the fit 
was poor. These results suggest a connection between 
the dominant power source in the system, and the mor- 
phology measured. 
Bussman et al. (2009, submitted) found similar results 



frame near-IR SED bumps (|Farrah et al.ll2008[ ). The 
morphologies are consistent with this picture. DOGs 
with power-law SEDs, thought to be AGN dominated, 
tend to show higher concentrations and small physical 
size. An AGN may be contributing to the centrally con- 
centrated light in power-law DOG morphologies. Of the 
five DOGs with the highest concentration (DOGs 2, 3, 
8, 10, and 11), GALFIT finds evidence for a significant 
point source (PSF/Sersic flux ratio greater than 18%) in 
four. DOGs with SED bumps, thought to be star for- 
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from HST ACS, WFPC2, and NICMOS imaging. Their 
paper showed that the luminous DOGs (/„(24) > 0.8 
mjy) have higher concentrations than other samples of 
z = 2 active galaxies including Lyman Break Galaxies 
and sub-mm galaxies. The high concentration measures 
were more pronounced in the NICMOS 77-band than in 
the optical (rest-UV) data. They posit that at longer 
wavelengths, the central AGN is able to contribute more 
light than at shorter wavelengths where dust obscura- 
tion is most pronounced. The AO X-band images of 
the luminous DOGs (/„(24) > 0.8 mjy) are even more 
compact than the HST i7-band (< Cao > = 3.4 vs. 
< Chst > = 3.0), further evidence of dust obscuration 
at shorter wavelengths, and higher AGN contribution at 
lo nger wavelengths. 

Das vra et al.l ( 2008) also show a morphology SED con- 
nection in a set of z = 2 ULIRGs (not necessarily DOGs). 
They find a trend of galaxy size vs. 7.7 micron PAH 
strength. Again PAH feature s are s ignposts of ongoing 
star formation. iDasvra et al.l (|2008l ) found that z — 2 
ULIRGs with strong PAH tend to have larger sizes than 
those with weak PAH. This is similar to our result that 
the bump DOGs have larger sizes than the power-law 
DOGs. 

4.2. Mergers as Triggers 

By analogy with local ULIRGs, the images may re- 
veal whether the DOGs are triggered by mergers. Mel08 
demonstrated the wide range of DOG morphologies. The 
3 systems discussed in Mel08 each exhibited different 
morphologies, an exponential disk, a de Vaucouleurs pro- 
file (suggestive of an elliptical galaxy), and an unresolved 
source. No obvious evidence of an ongoing merger (e.g. 
tidal tails, or double nuclei) was seen. 

The majority of DOGs in the current sample, 8 of 14 
that were fit by GALFIT, have resolved single Sersic pro- 
files with Sersic indices n < 2, suggestive of disk galaxies. 
In the Bussman et al. (2009, submitted) HST sample, 
the fraction of disk-like systems was even higher (27 of 
29). Only 4 of 14 DOGs in the AO sample have resolved 
profiles with indices, n > 2, suggestive of elliptical galax- 
ies, and one of those is better fit with a disk+PSF rather 
than an elliptical profile. Of the remaining three DOGs 
in the sample, two are consistent with being point source 
dominated, and one was too low S/N to fit with GAL- 
FIT. 

In contrast with Mel08, this sample contains two ex- 
amples of systems with multiple resolved components, 
circumstantial evidence for merging. DOGs 10 and 13 
show two resolved components with separations on the 
order of 1 kpc and total flux ratios on the order of one 
to five (possibly minor mergers) . They also both appear 
to contain an additional unresolved component, sugges- 
tive of an AGN. DOG 10 has a power-law SED and a 
very high concentration. DOG 13 has a bump-like SED 
and is very diffuse. These systems suggest that mergers 
may trigger both AGN and star formation dominated 
DOGs. However, the majority of DOGs in the sample 
do not show evidence for ongoing mergers (i.e. multiple 
resolved components). This conclusion was also reached 



in Busmann et al. (2009, submitted), which had only 5 
of 30 showing evidence for an ongoing merger. 

As Bussman et al. (2009, submitted) point out, the 
DOG merger fraction is much lower than in the local 
ULIRG sample which contains 35% with double nuclei 
at separations larger than 2.3 kpc (easily resolvable in 
both the HST and AO samples). The DOG merger frac- 
tion is also much lower than the merger f raction in the 
genera l z — 2 ULIRG population, in which [Dasvra et all 
(2008) found evidence for interactions in up 50%. Thus, 
if mergers are required to trigger DOGs, DOGs may pref- 
erentially exhibit the DOG photometric selection criteria 
only after the merger has occurred. Alternatively, the 
very red optical to IR color selection may just last longer 
than than the morphological evidence of the merger. It 
is also possible that mergers are unnecessary for the pro- 
duction of DOGs. 

4.3. Summary 

The morphologies of the DOGs presented in this paper 
appear to be correlated with their photometric proper- 
ties. As a result, a consistent picture is emerging for 
the nature of these extreme systems. The most lumi- 
nous, tend to show power-law SEDs and small, highly 
concentrated morphologies. Both of these trends can be 
explained by strong AGN activity. The less luminous 
systems tend to show SED bumps in the IRAC bands. 
They also tend to be larger and more diffuse than the 
brighter systems. The trends seen in the bump sources 
can be explained by a star formation dominated power 
source. Mergers may trigger both the AGN-dominated 
sources and the bump-dominated sources, although it is 
not clear if mergers are necessary to produce DOGs. 
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